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ABSTRACT 


Relationships relating power, magnetic field, current density, and 
ampere-turns in terms of certain dimensionless factors are summarized 
for many types of coil geometries and current distributions. A num¬ 
ber of plots of these factors are presented. 

The field homogeneity in magnet structures is presented in terms 
of a series expansion about the origin utilizing Legendre Polynomials. 

A number of tables to facilitate design of homogeneous fields are 
presented. A method of achieving homogeneity in long solenoid struc¬ 
tures by the use of determinants is discussed. Expressions for the 
axial field from uniform and radially varying current density coils are 
given. 
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PART I 

RELATIONSHIPS BETWEEN MAGNETIC FIELD, POWER, AMPERE - TURNS 

AND CURRENT DENSITY 


It has been customary la th8 literature on magnet design to write the relationships 
between fields, power, ampere-turns and current density in terms of a number of dimen¬ 
sionless factors. In this section we will treat in summary form twelve cases of coll 
construction (Cases I - XII). This work is in part a summary, and enlargement of parts 
of the following three papers: F. Bitter, R. S. L Vol, 7, 1936 Part II; W. F. Gauater, 
AIEE, Fall General Meeting, Chicago, October 1959; F. Gaume, Journal de Recherche 
du Centre Nationale de Recherche Scientiflque, No. 43, June 1958. (English translation, 
Lincoln Laboratory, M81-12, part I and H, Dr. H. H. Kolm.) 


A. The most commonly used dimensionless factor is the G factor first suggested 
by Fabry, which connects a given magnetic field with the power required to produce this 
field and is a function of the type of current distribution and the geometry of the coil. 


H =» G 



(i) 


with W - watts 

p = resistivity in ohm - cm 

Ql = Inner radius in cm 

X - fractional volume of conductor 

The space factor A is assumed constant throughout the volume. If A Is a func¬ 
tion of x and y, core complicated formulas are needed. 


B. A seco v dimensionless factor J has been suggested by W. F. Gauster and 
connects the current density at the innermost winding of a magnet to the power of the 
magnet. J is ah:. .. function of the current distribution and the geometry of the coil. 



C. tv. ....j ..vtlcular type of coil a maximum G factor can be found which will 
give the most magnetic field for the least power used. However, in many types of coil 
construction this maximizing of the G factor may lead to excessively high current densi¬ 
ties at the innermost winding which may make cooling the coil difficult. To illustrate 
this point a ratio, ~y , can be formed which compares the maximum current density in 
any coll to the current density In a uniform current density coil, where both have the 
dame Inner radius and both develop the same field. Quantities for *Jha uniform case are 
designated as primed: 
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and using (3) 



^ _ J g;^ 
jf = r G 



This ratio indicates the difficulty of cooling and supporting any coil compared with 
a coil of uniform current density at the same field and bore. 






D. The ratio in equation <3) is also useful and represents the power saving In 
using any current distribution compared with uniform current density when generating 
the same field in the same bore. We give it in the more useful form. 

W JQ '\2 

W' U j~ P (5) 

These dimensionless factors and ratios are presented In Table I for nine types of 
magnet construction and following Table I for three additional more complicated magnet 
constructions. 

One of the most Interesting new coil constructions presented is that of case V. 

This coll was first suggested by F. Gaume and achieves a relatively high G factor with¬ 
out excessive current densities and can be quite simply constructed of plates in the 
usual Bitter magnet style, with the introduction of axial current variations by mak in g the 
plates thicker towards the ends of the magnet. 

A number of plots of the geometry factor, G, against CT and 0 are given in 
Figures l through 6. Cockcroft's curves for case VIII are reprinted in Figure 1, Bitter's 
curves for case VI in Figure 2, new curves for cases VUl and VI in the small CT and 
0 region are given In Figures 3 and 4. New curves for case V (Gaume's axial current 
variation) and for case X (the ellipsoidal shell) are given in Figures 5 and 6. The 
variables CT and / 3 are defined: 
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E. It Is sometimes instructive to write the magnetic field in terms of ampere- 
turns rather than power. This irises from the practical standpoint that it is difficult 
to always predict the resistance of a magnet exactly, especially as a function of power, 
and one may wish to know how much field will be produced for each ampere available in 
the power supply. This relationship between ampere-turns and field can be related to 
the J factor and the G factor an! miscellaneous CT and /3 terms and numerical con¬ 
stants. The relationships are given here for three cases: Case VI, Case VII and Case 
VIII. All the others are easily derivable by Integrating the current density over the coil 
to obtain the total current in terms of J and W. Total current I and the factor J can now 
be substituted for W in equation (1), and the ampere-turns related to the field. 

Case VI: n f i *1 
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5*. Useful geometry factors can be written for very long solenoids (long enough 
that end effects can be neglected). 

For a constant current density coll the ampere-turns per unit length can be 


written: 






















NI ( (W/2b") A \ ,/J / 2(tr-J) V /s 

2b “ [aTT~) 

where W/2b is the power per unit length (watta/centimeter) 

For the disk construction, with a radial current distribution we write 



Rewriting (9) and (10) in terms of the field in gauss and power per length as watts 
per meter, we obtain , . 

* r F % 4 M mm ^ ^ 


where 


it / __\ (£ 7 rV fz f W (watts/meter) X Y /2 

H (gauss) -gj- G s /3 (ohm cm.) 7 

G 5 (uniform) = ^ 1 


G5 (radial) = ^J£m.< 3 *) 1</2 


These factors are plotted in Figure 7 and it is interesting to note that below an CC 
of 2. 5 there is negligible power saving with the disks, but at greater CC , the saving Is 
increasingly pronounced; at an CC of 10, 45% more power would be required by the uni¬ 
form wound solenoid for the same magnetic field. 

G. It Is useful to use the G factor notation when dealLg with split coils (coils 
separated by a gap). If we take a pair of coils with parameters as defLned in Figure 8, 
we can write the field In the gap by simple addition and subtraction of coils. 



Figure 8 


(ii) 










G; Gj (o*, ^3j) , ^2 = Ga (a*,/? 2 ) 


This relationship (11) holds for any current distribution having no axial variation. 














(Mtr 






































^*—*•"*• ~ «p»wt- ••wnn*w’<r*' 




!>l 


a 


a 




15 e 

r. 


!» 

* 

! 

t 


—Ji -‘' .-V — V.j 


CASE X (FIGURE 6) 

Coil Shape: Ellipsoids! shell, CC = constant 


cc - ap ya t 

yj „ ^ma* 





Geometry Factor: 
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sec 


(a z ~j3 z ) l/z 

Gmax =* .204 at O’= 6, ft = 2 

Current Density Factor J: 
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J/2 


^ ~ 1 4 -r/J in fl¬ 
at cr-6, >0*2, Gmax « 0.204 
and J mQZ -0.149 

Current Density Hallo ~y for same Field referred to case VIII.- - 1.31 
Vowel Hallo P for same Field referred to case VIII - 0. 763 
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CASE XI 

Coil Shape: Cylindrical Bore, tapered ends, non-coinciding apexes. 


slope k 



b =■ ky-t-o 
k = slope 


Current Distribution: 


Geometry Factor: 


t - io/ y 


i__ r k 

5 kCcr-O+y-kia* Lk 2 +1 


- ((k"+l , )a' 2 f2ky<r + y 2 ) , / 2 + ar (k z 4l) + y0^k z +l) 

((k 2 +1) + 2ky+y 2 )'/2 4 (k 2 +t)'/ 2 + y(^/ k z + A' /Z 

+ o m (cr ^k z ^l)42ky+y z ) l/z + y + k V| 

\ <Xk z fl) a: 2 -fSkycr+y 2 ) , / :2 + y ^ktz/J 
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G is a maximum when k = 1 and CC =4.5 which reduces case XI to case VII Current 


Density Factor J: 
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CASE xn 

Coll Shape: Cylindrical Bora, tapered ends, non-colncldlng apexes. 
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slope k 


Current Distribution: 


Geometry Factor; 


i- w h , _±qL_ 

/b ky+c 


ir'' z k 17 I \ 

5 (k*+l)^' L\ k (<r-l) + yi/n(X±2<_j/ * 


\p ( y? vA 
(0c*+ 0 a* 2 + 2kycr f y 2 ) +cr (1^+0 + y ^ 1^+?/ 

((k 2 +i) +2k^ fy 2 ;‘ /2 +(k 2 +l) !/2 +y 


G la a maximum when k = 1.0 and CJT =4.5 which reduces case XII to case VII Current 
Density Factor J: 

( 1 V /2 

J_ l 4Wkkr-tHy&^k+y 
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PART n 

HOMOGENEOUS MAGNETIC FIELDS 

A. 

A magnetic field can be written as a power series involving I^gendre polynomial# 
expanded about the origin of the field. If this origin is on the axis and on the plane of 
symmetry the expansion will have no odd terms and the axial component Hz and the 
radial component H p can be written: 

STzh)J K n ' 2 M ^ P 2n-2 («) 

Hp (r '^ J Tzhji [ h ^" 2)(z -°>] r2n ‘ 2 p 2„-2 w » 


U- cos & 


In these equations the coefficients are defined as follows: 

u (Zn-2), d 2n ' 2 H 2 (z,o) 

Hz (ZjO) = -- 


dz 


2n-2 


Z = 0 


and P2n-2( u ) are the Legendre polynomials tabulated in Table 2. 


The importance of this approach, due to W. Garrett 1 , is that the field at any point 
on the axis of a system of coaxial coils can be given explicitly In simple algebraic form. 
From this, the derivatives of equations (1) and (2) can be obtained, and so a complete 
expression for the field in the vicinity of the origin can be obtained. If sufficient deriv¬ 
ative terms are used the field can be found anywhere within a sphere whose center lies 

' i 

at the origin and which does not extend far enough to iDcludo any corners of the coil. How¬ 
ever, out to within a few percent of the Inner radius of the coil, the field can be found 
quite accurately with only a few derivative terms. Compensation can be designed for each 
component of the field separately. In general, if we have a coll system of the assumed 
symmetry' having N variable parameters such as size, or ampere-turns, it is possible to 
specify N coefficients in the expansions 1 and 2. 


For example, to examine the standard Helmholtz pair of coils, where the coils aro 
so spaced that there Is no second derivative of the field at the origin, we can think of 
the coil system In the following way: 


1 . 


M. W. Garrett, J. Appl. Phys. 22, 9, Sept. 1951 













I 




(a) (b) (c) 

OL b “G 

Figure 1 

As shown in Fig. 1, we consider a single coil having length equal to the coils and the 
gap separation, from which we subtract a second coil of the length of the gap. The gap 
coil Is thought of as having a current in the reverse direction to the original long coll. 
We now merely need select the length of the gap coll so that it will have a second deriv¬ 
ative equal and opposite to that of the original long coll and the two will thus cancel at 
the origin. 

There are many other combinations which would also cancel the second derivative 
at the origin, for Instance, the hypothetical gap coil with its reverse current could have 
been somewhat longer and had less reverse current in it, or shorter and had more re¬ 
verse current in it. Roth of these solutions, of course, lead to a continuous solenoid 
with no gap but with an altered current density in the center section. These various 
solutions will have both an effect on the amount of field in the center of the magnet and 
on the size of the higher derivatives. The greater the distance over whdch one compen¬ 
sates, the more parallel will be the flux lines from the compensating coil, hence the 
flatter will be the flux distribution at the origin, and hence the smaller will be the sum of 

the higher derivatives. This compensation may not be compatible with obtaining the 

2 

highest central field for the least power, however. 

B. 

Equations (1) and (2) can be rewritten In the following form.-; 

H* (r,o) - H z (o,o)[n-£ z £«) + <£*(jj;) P 4 («) + ~] » 
H p (no) - Hz (oio) Jo+£ z (ai) Zp 4 ( u ) + £a (£) P 4( u ) + "'j (4) 



2- R. S. Ingarden and J. Mlchalezyk, Bulletin De L’Academie Polonaise Des Sciences 
(Serle des eel. math, astr., and phys) VIII, 5, 1960. 
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The expansion has been normalized to the Inside radius of the coils a^. A tabulation of 
the £it coefficients for three useful geometries is given in Table 1. The geometries 
treated are (i) a cylindrical thin current sheet, (li) a finite solenoid with a uniform 
current density, and (in) a finite solenoid with a radial current distribution (i^» %>• 

The value of the first eight Legendre polynomials and their derivatives are given for con¬ 
venience in Table 2. From equations (3) and (4) the magnetic field can be found at any 
point near the origin of a magnet, at a distance f from the origin and at an angle O - 
The derivative terms In equation (3) have their largest value for O = 0 (P (1) =1 for all 
n when O’ ~ * and cos O = 1)- The maximum value of the terms therefore can be ob¬ 
tained from the simpler form of (3) where O ~ 0 smd Z =• V 1 

H z (2 ( o) =* H z (0j0) jj 4- 64. 4- — J (5) 

c. 

The most commonly used method of achieving some degree of homogeneity In coils 
is to use the spaced helmholtz pair. Garrett has compiled a number of tables dealing 
with this method of compensation and since these tables have not been published we are 
making them part of this report- These tables are for constant current density only. 

The first table. Table 3 is used to determine the proper spacing between finite coils to 
achieve cancellation of the second derivative at the origin. The second table, Table 4 
is used to determine haw far from the origin one can move before the field deviates by 
more than 0.1 percent of the field at the origin. The third table, Table 5 is particularly 
useful and can show how accurately the coils must be spaced in order to cancel the 
second derivative to the desired degree. The two remaining tables. Table G and Table 7 
give the resultant fourth ard the sixth £ coefficients x'emaining after the proper spacing 
of the coils. The tables are used with the notation of equations (3) and (4) except that 
Garrett has normalized tLe expansion to the mean radius of the coil rather than to the 
inner radius as In Section 13 above ( (~n) rather than (j?L ) where ) ). 

To illustrate the use of Garrett’s tables, we work out the following examples: 


Assume we have two coils, each with an CT - 3 and a ^ 1 and we wish to sepa¬ 

rate them to cancel the s-econd order derivative at the origin. 


X - 

A =* 


4/3 = 

a-i ] 

2(g-i) 

cr b l 


1.0 

« 1.0 


From Table 3, 


K - -s*2ro 


AX 


A X 

atO + crT 
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:y 

fiafoVi. 


therefore 


- 2.35040 - 2/3 4 -^ 

ai ' ai 


. ;'3040 


=* . 181 and ~ =. 724. 

ai 


This means that 


and the spacing should be 

M 

at 

Using Table 4 we see that V - ~~ =. 181 or qi (i-hcr) 
out to a distance along the axis of 72.4% of the inner radius, the field is within 0.1% of 
the -field at the origin. Using Table 5, we can examine the effect of improper spacing. 

H= H 0 (i-i-Q (4) -(^) 2 + £4 ( 0 + -■) 

where Q is the value in table 5 and d/a Q is the displacement error In locating one coll. 
For the above example, we can find the error that will make the uncancelled 2nd deriva- 
tlve error equal the fourth at ( — ) = 0.1 

Q - 4.21 (at X-1,A = t) 

£ 4 = -922 from table 6 
4„2t ^ x 10" 2 * -922 x 10 ' 4 - 

A = 2.18 x 10- 3 

a 0 


D. 


When maximum homogeneity and not maximum efficiency is of the greatest impor¬ 
tance, the higher derivative's must he handled. Since the number of derivatives that can 
be cancelled depends upon the number of variables, we must now pick solutions which 
have more variables at our disposal. For Instance, If we wish to cancel both the second 
and the fourth derivative of the field we can uuperirnpo.se two colls on each other each 
with two geometry variables or each with oue geometry variable and one current variable. 
To proceed, then,we would pick the length of one coil and the current of one coil and the 
variables then v/ould bo the length of the second coil and the ratio of the two current 
densities. Garrett has calculated a set of so called sixth-order solenoids (no second or 
fourth derivatives) by maxing the following choice of variables: lie constructs an over¬ 
wound end solenoid by winding N layers of turns over a certain central region of the sol¬ 
enoid and 2N layers over a certain section of both ends (l.e. a notch on the center 
section of the o. d. of the coil). Ho now solves uniquely for the length of the overall 
coll ami the longth of the notch. The results of his calculation are shown in Table 8 for 
several thicknesses of coil. 

There are oilier choices of two variables that one cun make to construct sixth-order 
solenoids and Table 9 illustrates several methods. Table 9 was constructed for the 
specific example of a ra'her long solenoid. Several general remarks can be drawn from 
the table. Probably the most Impeltant of which is that since U 10 solenoid is very long 


e 


• • 












the job of compensating it for high homogeneities is much simplified, the higher derivative 
having been quite small to begin with. We notice, however, that in all but cne case In 
the table, cancelling the second and fourth derivatives of the field, have increased the 
sixth derivative. The single exception to the current sheet of the o. d. The current sheet 
on the o.d., however, requires more power than any of the other methods of solution. 
Perhaps the meet attractive method of solution is that represented by the last column, 
that of current sheets on either end of the i.d. of the coil. This method takes little 
power and leaves the center section of the magnet open for further compensation if desired. 

Cancellation of higher derivatives than the fourth is of course possible, but is often 
of decreasing practical Importance because of the unrealistic accuracy necessary in the 
cancellation of lower derivatives, if the leftover uncompensated part of say the second, la 
to be smaller than the sixth derivative or eighth derivative, etc. It probably makes sense 
to design colls to cancel the second and the fourth derivative and use a small auxiliary 
set of coils driven from the separate power source to tune out any random errors which 
creep in. Some of these random errors are of course going to be nonsymmetrlcal and 
will give rise to odd derivative terms. It la quite possible to construct a set of auxiliary 
compensating coils which will operate independently of each other. One coil can be 
designed to have a second derivative and no fourth and a second coll to have a fourth 
derivative and no second, etc. Each one can then be tuned independently of the other. 

E. 

We have constructed a table of £. coefficients up to the eighth order for coils of 
varying <X and J3 with both constant current density and with current density equal 
i - . The tables also Include the G factor and the multiplication of the G factor 

and the & coefficients. The tables are to be used with the expansion of equation (3) 
and (4) where z Is referred to the inner radius Qj ; & Is tabulated as D for uniform 
currents and a 3 A for the radial current distribution. Table 10A gives the coefficients 
for the uniform current density case for several CT 3 and for up to 10. Table 10B 
gives the coefficients over the same range of CC and /3 for the radial current distribu¬ 
tion case. Tables 11 A, B, C, D show more detail for the constant current and the 
r adlal current cases in the small range. 

The use of the tables is probably best illustrated by the following example: Let 
us assume wo have a coll we wish to separate for puiposcs of access and we want to 
know the resultant Held In the gap and the resultant homogeneity (see Fig. 2) 
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If we now subtract from the field expansions for a coll of length 4b+jt^, a coll of length 

-t „ we can write: 

£ 

H gap - (Ho t “Hog) + (H2t“H2 g )^)F(H4 t -H 4fl )f^+- — 

-■(sHs 

i . 


A Y /2 f7„. .L „ eo^ ,/a . .V*’ 


II-W 


W* / ' l GtE2 W~ £z ’ 


2 

+ 


V t4 ■^ t G/~‘ £ 4)^J + — 

where the primed quantities refer to the error coefficients for the gap. TaViag a spe¬ 
cific geometry, let the current he corstant (5 = D). 
let 

fit “ 4 

a* =» 3 

/?g = 0.2 

What is the field in the gap and the second error coefficient, and what would be the 
second error coefficient if there were no access gap? 


•> 














The D? listed for tT = 3, /3 -4 in Table 10A is - .1415 x 10*, therefore the second 

c* - 2 

error coefficient of the whole coil would be £2 = -1.415 x 10 where 


H - H 0 (I + £ z + "-) 


To solve for the split coil, we use (5) and the Gj . Dn columns in Table 10A 

•V* OTIH -0.08639( &&)*) + 

+ 0.2235 * 1(T> 0.5004 x 10^ "*J 

H 9°P ( f 03) [• 1587 + °' 8965 * 10 ' 2 (li) 2 ] 


this gives a central field of 


Hgap (o,o) - ( 0 J 4 3) 

and If the second order coefficient Is written In the following form: 

H=Ho (1 + £2 (f,)%-•) 

the coefficient for the above example would be 

€ z a = 6 46 x 1o -Z 


.1687 


The loss in field for this case has been 9. 5% and the second order error coefficient has 
been Increased by a factor of 4.55. Lf we had been trying to homogenize the field by 
separating the colls, we note that we would have separated them by too much. To ”ie 
these tables for synthesis (l. e. cancellation of derivatives) It would be necessary to use 
trial and error. In the manner of the above example. 


It Is Interesting to note the oscillatory behavior of the £ coefficients, particularly 
Lu the small yQ region. The second derivative Is always negative, the fourth derivative 
goes through zero once, the sixth derivative goes through zero twice and the eighth 
derivative three times. This means that there are colls which can be built with no fourth, 
no sixth and no eighth derivatives. This may be useful in the design of Independent com¬ 
pensating colls. 



















..'Mt Mmh> u. ff— ... *2Li&i*. Anir' . 




F. 

If one wishes to compensate a very long solenoid over a considerable length of the 
axis such as a multicoil structure used for plasma experiments, expansion of the deriva¬ 
tives around the origin has limited usefulness. One might wish to use the following pro¬ 
cedure instead; namely specifying that the field at any number of Intervals along the axis 
be constant and then solving for ti e current distribution in the coils that would make the 
field be equal at the specified points. In the case of a multicoil structure one could ask 
that the total field in the center of each of the coils be equal to a constant and then, solve 
for the current density required in each of the coils to make this true. The same proce¬ 
dure could be used with a single long solenoid where the solenoid could be broken up into 
an arbitrary number of pieces. 


The most useful method of solution of this type of problem is basically the following: 
taking any single ceil we first find the field at the center of that coil and at distances 
away from the coil represented by the center lines of each of the other coils, for a given 
current. Then by properly adding up the colls in matrix form and demanding that the 
field matrix times the current matrix be a constant we can invert the current matrix and 
solve for the individual currents required in each of the coils to achieve the desired re¬ 
sults. The method of solution and a useful notation are developed in the following example. 
Mr. R. Bradshaw of Conesco, Inc., Arlington , brought this method of solution to our 
attention. Consider the following system of coils arranged so that they have a common 
axis. They may in general be of various sizes but for illustration consider them to be 
all identical with the same separation. For ease in calculation choose an even number 
of coils so that we have a plane of symmetrv passing between coIIb. 



in each coil 


Our object is to find what ratio of currents is needed in the various coils so that a uni¬ 
form magnetic fie< • ’ is obtained throughout the Interior. If the coils are thin and close 
together, to a approximation uniformity is obtained by requiring that the field at 

the center of ear. coil be the same. The field falls off rapidly on each side of the 
center of a . ■' < > and we may construct a matrix the elements of which represent 
the contributions o. i'v colls at their respective centers. Let Fij be the field at the 
center of coll i duo to the field of coil j . For example, for 4 coils, 2 coils on each 









/..y. ? ;. ... ; \ . . •; V ^ 


■*• “* »« • -nr ■»> y.i-^y-ff«i'"“^TWV , y*« »-■%*»-*• 


• • 


side of the plane of symmetry, we have 

I 

F 1S p i2 | F|3 Fj4 

Fzt F 22 jF 23 F 24 

^31 *32 3 F 33 F34 

I 

F 41 F 42 I F 43 F |4 

I 

These matrix elements may be expressed In % of an Infinite solenoids! sheet of the same 
mean radius, or whatever one desires, since we need only ratios. 

Fij Ij - ?> where *J\ is the field at the center of the i *^ 1 coil or 

: 1 ms si 


• 9 



If 


7, 


A 2 


7 2 


I 3 


7 3 


I 4 


74 


F 31 F 32 I F 33 F 34 I 3 y 3 

Ft2 | f 43 F44 I4 y A 

If one now specifies ~*ie the Fij matrix can be inverted and the Ij matrix 

found. 

If the current in the coils are arranged so that If~l 4 and I# "I 3 then It is convenient 
to partition the matrix as follows: 


f _ hi Fiz f _ F,j F ,4 F 3i F 32 F Fm 

*" ~ IfaF« ' f,z_ F23F24 ’ F 4 ,F 42 • ^ f 43 E)4 


JFmFm 


• • 


• m 


• • 


• • 


• • 


: , 13 . 

I 4 * 




__ 


• • 


fn ftt 

ti 

- 

and then x 2. - 

^21 ^22 

*2 


where k =*J 


01 1 , 


fn 

fl2 

U 

- 7{ 

^21 

?2Z 

ki s . 

k?j 


• • 


• e 


it follows that 


(fn + fjek)!, = y\ 

















rv 


-Ji 




lukij .'Jl 


_»aXk.4 A. jjj.UJj fr W^A 


■■'. V A.fc.»uA I l£/iz£~ 


Since we desire the fields at the centers of each coil to be a constant 
J'X ~ 3*2. ~ constant = C whore C 3 C | j j 

*.e., (fu + = C 

It Is clear that if -(f^t + hz^c and this gives us the desired ratio of currents. 


As an example consider the field on the axis due to a system of four current loops 
of radius y a 1 Cm 



li- U 
I 2 "l3 


Figure 4 


Separation between centers of the coils is 1 cm. We need calculate the field of one 
loop at the center and 1, 2, and 3 cm. on the axis away from the center. For a single 
loop 


then 


H x 


2 rj y 2 

10 (x 2 Hy2)5/2 



0.628 0.222 

0.222 0.628 



fick 


0.056 0.019 
0.222 0.056 


|0 1 
|l 0 


X 

Hx/i 

0 

.628 

1 

.222 

2 

.0562 

3 

.0199 


0,056 0.019 

0.222 0.056' 


0.019 0.056 
0.056 0.222 


fii + ffzk. 


_. A 

0.628 0.22Z 


0.019 0.056 


0.647 

0.278 


0.222 0.628 


0.056 0.222 


0.270 

0.850 


We now need the inverse of ( fll + flZ^O 


By definition, tha Inverse of a matrix A with matrix elements djj Is obtained if for 

each element (Q|j) ^ — Li.?— l ^ 2L J Aji is the minor of the determinant A. 

det A 
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if A — I a,f a i*| then A -1 - 

I «2I a 22t 


For 2x2 matrices this is very simple: 


( a n) f (a 12) * 

( a 2l) 1 (^Z2) * 


A -1 - 


det 


|_ / a 2Z -0-12^ 

tA \~&21 Q,i/ 


det (fij + fizk) - (.850)(.647)-(278) 2 =. .473 

(fw+W-i 


.850 -.278 


-.276 .647 



1.797-.587 
-567 1.367 


! 

1 

If 

1.797 

-587 


1.210 

• • 

{ 

h 

" -.587 

1.367 

| C 5=3 C 

.780 


•*. If » c(l210 ) I 2 - C (0.7Q0) 


The current In coil 1 must then be 1.55 times greater than the current in coil 2. 

Ingarden and Michalczyk^ have also considered the problem of homogeneities over a 
certain length of the axis rather than strictly at the origin. They have used the require¬ 
ment that the mean square deviation over a certain interval be minimized. Their notation 
is useful, but solution of the equations, which are transcendental in character, presents 
a number of problems. 

G. 

It Is often of Interest to know the field along the axis of a solenoid both inside the 
solenoid and outside the solenoid. The following two equations can be used to find the 
field along the axis of a uniform current solenoid and a radially distributed current sole¬ 
noid (Bitter discs) respectively. 

(7) 

kZM/2 


(I) Uniform current lc« %, 


,, TriXaiffl (k-/3)4-( r ct 2 +(k-/3) z ) l/ ^ 

H " 5 L <lc+/»+(«*+(k +fi) 2 ) l/z 


Q< /3) + (lb(k-/3)^)^ 1 

(k+£H04-(k+/3; z ; v2 J 


(ii) Radial current hjj. 


k= £- 

K ai 


iruXaX,. ^ 0 I cc+CgZ+Ck+fl) 2 )'* 

H “ ~5-!(k+/3)fin|-J f(1t ( k+/ ,j2)i/2- 

arHar + Qt - epF * } 
'M** i +(t+(k-0)V' a J 


( 8 ) 
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Tabulated values of equation (7) for a wide variety of coll shapes have been published by 
D. E. Mapolher and James Snyder in a University of Illinois circular #66 entitled. "The 
Axial Variation of the Magnetic Field in Solenoids of Finite Thickness". There in no 
available tabulation of equation (8) 

Equations (7) and (8) can be written in a much simpler form utilizing the following 
angle notation: 



02 

a, 

b. 

at 




- a 


/3 


- ~y 

CL| 


i) uniform current density 

TriiXai 

H =» —_— 




Hi+W ' H^¥) 


(9) 


II) radial current density p 

l -ln 


11 _ TitXa t 

H--g— 


tan ^ tan ^ 


tan • 


#ii 


fan 


( 10 ) 


2 ™ 2 

General expressions for the off axis fields arising from a loop of current can bo 
written In terms of elliptic Integrals of the first and second klnd.^ They can be written 
In the following way, where the parameters are defined in Figure 6. 


H 


Ik [r(kH ( ^-^- 22) 


Ql/ ,| r(k) ^ -r - 7 sW 


h = -^- 

Ql/2 


T * (a z +p :£ + z 2 ) . » 

~F(k) + v ^- >- 3 (k) 


'] 


(ii) 


( 12 ) 


3. 


W. R. Smythe, STATIC AND DYNAMIC ELECTRICITY, p. 266, McGraw-Hill, New 
York, 1950 
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Q =■ ( a+/o) z +z z 

1 _( 4 ap Y 2 

*"Wj*+W 


F(k) 

s(k) 


Elliptic integral of 
the first kind 


E(k) 

1-k* 


where E(k) - 
elliptic integral 
of the second 
kind 
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Using equations (11) and (12) the field at any point produced by a finite coil can be found 
by dividing up the coil Into loops and summing all the contributions. A number of labora¬ 
tories, including our own have written simple computer programs utilizing equations (11) 
and (12) to find the fields in and outside the conductor volumn. 


On-axis and off-axis fields from any solenoid with any radial current distribution 
can be calculated by means of a lengthy set of tables published by the Oak Ridge National 
Laboratory Report #ORNL 2828, entitled, "Tables for a Semi-Infinite Circular Current 
Sheet, " by N. B, Alexander and A. C. Downing. These tables are for the field around 
the end of a semi-infinite current sheet and the coil in question is built up out of these 
current sheets In the method shown In the introduction of the table. 


Attention is also called to a tabulation by E. E. Callaghan and S. H. Maslen, NASA 
D4G5 "The Magnetic Field of a Finite Solenoid". This paper plots graphs of the radial 
ami axial components in and around current sheets of various lengths. While not as 
accurate as the Oak Ridge tables, it is very useful for approximate results. 


A final paper of particular relevance is that of G. R. North, "Some Parameters of 
Lumped Solenoids", (Oak Ridge ORNL-2975). He derives an expression for the field 
ripple resulting from the separation of coils in a long multicoil structure (end effects 
neglected). The equation is as follows (see Figure 7) 




Sin 


2b 


2brr 


• COS 


2ttz 


2 i*" 1 
s’ 
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1/2 - 

or e 


2-needs 
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(13) 
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TABLE 1 

TABLE OF ERROR COEFFICIENTS 

THE £ COEFFICIENTS ARE CONSTRUCTED FROM THE 
FOLLOWING VARIABLES: 


c . = 7T2T a ~ “ s '" z 3i 


c? - = t 5 2 ' = cos2 ^ 



2 Ta, 



^4= 


cr 

si n 2 # n 


a z +/3 z 

10 

a 2 

— - = a: 

a i 

/3 Z 

cos 2 

£ = /3 

a i 

<3C Z +fi Z 

10 


« / cr4-(/3 2 +cr 2 ) ,/2 

M 1+(/? 2 +1) ,/2 


tan (% + ° 10 /2) 
tan(% +^"/2 ) 


38 im. [ or 


/?+(/3 2 -H) 

,0 -*-(>3 2 +cr 


1/2 \ 

5*7 


tan ( ^ 
tan ( ^ 


,0 /2) "j 

”/2; J 
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THE COEFFICIENTS ARE USED IN THE FOLLOWING SENSE: 
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H Z = H 0 (l + 4(j) 


P* W +<5 4(57j^(u) + 


• • 


• • 


u a-..\ • 


U“ COS# 
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TABLE 1 (continued) 


CASE 1: CURRENT SHEET (I = AMP TURNS/cm.) 


H (r } &) - cos# 11 4-^ 


cf 


£4 ~ ~ *A c f — 7c j 4 ) 


£, 6=-g i (^ c 2 2 ' 3,5 ^ +:! i § ) 


• • 


(a) P' ? ^ + (o') p4 ^ U ) + "* 


• • 
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CASE 2; UNIFORM CURRENT DENSITY SOLENOID 

HM“G(^)' / 2 =Ho 

H (r, 0 -) - H 0 ^1 + £ z (u) + £4 % ( u ) + 


1 („ 3/2 ^ 5/2 \ 

2/3 2 c 5 ( Ci ’ 3 ) 

= llvk s [ 0 ' I ° z+ f C 2 2 )-°" /2 (' + 1 C 4+ f Cz ) 
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TABLE I 


CASE 2 (continued) 


, 15 2 35 3 315 4\ 

, rr^s c 2-t e 2+T^j 

cr(i + i-c 4+ f c ,_ fc 3 + ^ c ^ 

' /l! ( 1 + t Caf + ff C 2 ^ 

^?-c 2 5 + ^c 2 6 j 

• 3/2 Af -i. -3^ . 15^ 2 , 35, 3 , 315 4 

-3 p+2 C ^+^8' C 4 + 76 C 4+li-C 4 

2079 „ 5 , 3003 .. 6 \~1 
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TABLE 1 (continued) 

CASE 3: RADIAL CURRENT DISTRIBUTION (**'%) 
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TABLE 2 

TABLE OF LEGENDRE POLYNOMIALS 
AND THEIR DERIVATIVES (EVEN NO’S) 


Po ( u ) = 1 u- co& ^ 

Po 00= 0 

/°z (u) = ~(3u a -l) 


:• 




/4 ( U )=-|-(6u) 

/°4. ( u )=-g-(55 u 4 -30u z +3) 

/°X (u) = -^-(l40u 3 ~ 60u) 

( u )“ )6 ( 23 ,u6 -3!5u 4 + 105u z -5) 
pi (u)=~(13S6u 5 -1260u 3 + 21C>u) 

Pa ( u )=j|g (6.435u« - 12,012u 6 f- 6,930 u 4 - |,260u 2 +35) 


p'a 00 = 755 (51.4S0u 7 -72,072u s + Z7.72Ou 3 -Z,52O) 
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AXIAL ERROR LIMITS 
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